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Chapter 1 
Introduction 



1.1 What are we made of? 

Where we come from, where we are going to and what we are made of are 
questions that have been asked by humans at all times. In former days, the 
only answers that could be given were speculations by philosophers or claims 
by priests. Even nowadays, we cannot answer any of these three questions 
in a satisfiable way. But we have shifted the emphasis of them a Uttle. 

The quest for our origin is nowadays shifted to evolution and its means, 
and / or finally to the big bang that astrophysicists beheve has happened 
about 13.7 bilhon years ago. No-one can answer the question where that big 
bang came from. 

Defining our destiny is something that imphes the need for forecasts; as 
our understanding for the other two questions grows, we gain more and more 
knowledge on how to forecast things. Still, a precise forecast on the world's 
state tomorrow will never be possible. 

The third question remains: What are we made of? We can neither 
give a final answer here. But we have come an enormously long way from 
the beheves of the oíd greeks, that everything consists of water (Thales) or 
maybe of the four elements fire, air, water and earth. In the 19th century, 
chemists found parts of the matter they considered indivisible, and they 
classified about 50 so called chemical elements that are built up from those 
indivisibles — the atoms. Nowadays we know 91 natural chemical elements 
and have synthesized several more; their number is about 114. We also found 
out that there exist differences between the atoms of a single element; their 
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masses vary. Atoms of different mass belonging to the same chemical element 
are called isotopes. We know about 2.500 different isotopes. 

The sheer variety of classes of atoms (and the fact that we can synthesize 
more) indicates that these are not as indivisible as chemists thought in the 
first place. Indeed, already in 1909 Lord Earnest Rutherford of Nelson found 
that atoms — he experimented with gold — consist of a very heavy, positively 
charged core, called the nucleus, and a very light shell, which is negatively 
charged. It was soon identified that the shell consisted of the same particles 
that make up for electricity, the electrons. The nuclei, on the other hand, 
have been found to consist of protons and neutrons, the so-called nucleons. 
Roughly spoken, the number of protons in a nucleus determines the chemical 
element, whereas the number of neutrons distinguishes different isotopes. By 
those discoveries the variety of 2.500 different indivisible particles has been 
reduced to 30- 

When trying to study electrons, protons and neutrons, further particles 
have been discovered. Besides particles that resemble light — the so-called 
photons — and similar particles, the new particles can be put into two classes: 
Those that interact with nuclear matter as nuclear matter itself does, i.e. with 
the so-called strong forcé, and those that don't. The first are called hadrons, 
the latter are lepton^. In the standard model of partióle physics, which is the 
state-of-the-art-theory of our knowledge about matter, the six known types of 
leptons (called electrón, muon, tauon, electron-neutrino, muon-neutrino and 
tau-neutrino) remain fundamental parts of matter. No experimental results 
indícate that they have an inner structure. 

This is not the case for the hadrons. Hundreds of hadrons have been 
directly or indirectly detected, and it turned out that they are not funda- 
mental. Instead, new degrees of freedom have been postulated and experi- 
mentally observed, so-called quarks and gluons. In accordance to the number 
of leptons, we know six quarks, called up, down, charm, strange, top and hot- 
tom, from which only up and down form neutrons and protons. Gluons are 
to the strong interaction what photons are to electrodynamics, the theory 
that covers electricity, magnetism and light. 

Quarks and gluons leptons and photons, considered fundamen- 

tal parts of matter. Their observation, though, is not very easy. Due to 

^For complction, it should be said that the electrons make up for the chemical properties 
of an element. but their number is dctermined by the number of protons in the nucleus. 

^The terms have greek origin and mean "strong particles" (hadrons) and "light parti- 
cles" (leptons), referring to the strong forcé and low masses, respectively. 
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the structure of the strong forcé and the underlying theory, the Quantum 
Chromo Dynamics (QCD), they can never be seen alone, but only in pairs 
(one quark and an antiquark) or in triplets (three quarks or three antiquarks), 
or any combination thereof. This phenomenon is called confinement. When 
scatterings between quarks happen at very high center of mass-energies, this 
confinement ceases to exist. This is called asymptotic freedom. It is sub- 
ject to experimental research to créate deconfined matter, i.e. matter that is 
asymptotically free. Since this matter consists of free quarks and gluons, it 
is called Quark Gluon Plasma (QGP). 

1.2 Studying Quark Matter 

The purpose of the largest experiments on earth, such as the Relativistic 
Heavy Ion CoUider (RHIC) at Brookhaven National Laboratory (BNL), the 
Tevatron at Fermilab and — soon — the Large Hadron Collider (LHC) of 
the European Council for Nuclear Research (CERN), is to probé the smallest 
parts of matter we know. At these facilities, protons (Tevatron and LHC), 
gold- (RHIC) and lead nuclei (LHC) are accelerated to velocities very cióse 
to the speed of light and coUided with each other. In proton-proton-collisions 
as the smallest possible system of coUiding stable hadrons, the fundamental 
forces and particles can be studied in a very clean way, whereas in the bigger 
systems (gold on gold and lead on lead) collective properties of the matter 
can be studied. 

The latter is the kind of physics that is addressed in this thesis. In systems 
with 396 or 416 initial nucleons (gold and lead, respectively) , it can e.g. be 
studied how the single particles react in connection with and surrounded by 
many other particles. Quantum Chromo Dynamics (QCD), the gauge theory 
for strong interactions, predicts free quarks and gluons in the infinite high 
temperature limit. It is hoped that this "deconfinement" can be reached in 
collider experiments. Indeed it is claimed to have been seen in experiments 
at the Super Protón Synchrotron (SPS) at CERN [1] and at RHIC. 

When studying heavy ion reactions, a major field of interest is the so- 
called phase diagram of QCD. Like in everyday physics, such a diagram 
shows in what state the matter is at given conditions. Unlike in everyday 
physics, the conditions are not given in pressure and temperature, but in 
(baryo-)chemical potential /ig and temperature T (see figure [TTT]) . Depending 
on the center of mass-energy y/JÑÑ of the collision and on the system size 
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T 



170 
MeV 




310 MeV 



Figure 1.1: The phase diagram of QCD. It shows the different expected 
phases. At /i = 310 MeV and T = O lies ground state matter. Note that the 
position of phase transitions is not known exactly. From [3] 



(determined by the nuclei involved) one can probé different regimes of the 
diagram. With very high energies and big nuclei the matter is heated very 
much (high T), but has few baryons (low /i^). Such matter is created at 
RHIC and LHC. Here, a second order phase transition is beheved to occur 
at a temperature Te ^ 170 MeV [2]. 

This phase-transition leads from a gas of normal, confined hadrons, the 
Hadron Gas (HG), to deconfined Quarks and Gluons, the Quark-Gluon- 
Plasma (QGP). This phase seems to behave like a liquid, according to claims 
made by the experiments at RHIC |1] . This is surprising, because it has pre- 
viously been thought of clS el perfect gas. 

The theoretical problem one has when dealing with nuclear matter is that 
QCD cannot be solved perturbatively at small energy-momentum transfers 
Q^. This is due to the at low large coupling constant. Unlike the elec- 
tromagnetic coupling the strong coupling varies a lot and actually falls with 
increasing Q^, which is the reason for the theory being "asymptotically free" . 
For lower temperatures (and therefore lower Q^), nuclear systems cannot be 
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calculated from first principies (namely QCD), but with effective models. 
Some of these are described briefly in chapter [2l Also, the concept of tem- 
perature itself introduces the need for many particles, and there is no kind 
of interaction in which many-particle systems can be described exactly. This 
is another reason for the need of effective models, even above the threshold 
for deconfinement. 

Of the existing models for dynamical descriptions of heavy-ion reactions 
(see chapter [2j), hydrodynamics is very popular. In this approach, the matter 
is described as a fluid. Henee, coUective effects can be studied easily. For ex- 
ample, collective flow is an observable intimately connected to hydrodynamic 
behaviour (for more detall, see section [23]). 



1.3 Supersonic acoustic sources 

In all fluids, there exists a mode for propagating weak, linear perturbations. 
Those perturbations are exactly what excites the human eardrums and causes 
US to hear — sound. The speed with which these perturbations move is the 
largest speed by that any mechanical stimulus can travel through a given 
body (which does not necessarily have to be a liquid) and depends on the 
material the body consists of. It might also depend on the wavelength of the 
perturbation, this phenomenon is known as dissipation and does not exist in 
perfect liquids. 

The speed of sound cs can be calculated from the Equation of State (EoS, 
see section 12.41 page [T2]) to be the partial derivative of the pressure p with 
respect to the energy density e 

4^1 . (1.1) 

When an acoustic source moves through a médium, the audible sound 
changes. A resting observer in front of the source will hear a higher fre- 
quency than is actually emitted, if the source passed her, she will hear a 
lower frequency. This phenomenon is known as the Doppler-ejfect. It can be 
easily understood if one considers the sound waves emitted by the source as 
being compressed in forward direction and elongated in backward direction. 

If the speed of the source is the same as the sound velocity, all sound waves 
ever emitted by it will reach the observer at the same instant — together 
with the source itself. One cannot hear no frequency anymore, but only a 
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Motion of front 



Figure 1.2: A schematic view on how a mach cone develops (in two dimen- 
sions): The elementary sound waves all add up at one line. The angle a is 
given by eos (a) = cs/v, where v > cs is the velocity of the source. 

(super)sonic boom. In the supersonic regime, when the source is moving faster 
than the sound it emits, the sound arrives after the source has passed. The 
waves form a cone — the mach cone, named after Ernst Mach, an Austrian 
physicist who hved from 1838 until 1916. The cone develops because there 
is a straight surface perpendicular to every elementary sound wave emitted 
at any instant, see figure [L2l 

Unlike normal sound perturbations, the amplitude of the mach cone does 
not decrease with the square of the distance to the exciter, since it is not a 
point, but only with the distance to the first power, since the exciter forms 
a straight line (this is analogous to the electrical field of a pointlike charge 
and a infinitely long, homogeneously charged wire). 

Measuring the angle of a mach cone will give insight on the speed of 
sound or the speed of the particle, if the respectively other is known. There- 
fore it might help to falsify an assumed Equation of State for the médium 
considered. 

Whereas in everyday experience, the velocity of the médium, say, air, is 
small in comparison to the speed of sound and its fiuctuations therefore are 
small as well, this is not the case in general. In nuclear matter, coUective 
motion may be as fast as and even faster than the speed of sound and vary a 
lot along the trajectory of the sound source. To study how much this affects 
the resulting shape of a mach cone is the goal of this thesis. 

In our special case, we consider a high energy jet, as might result in 
a partonic collision in the very early stages of a heavy ion reaction, and 
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examine the sound waves it emits and their way through a realistic médium 
that moves with a speed comparable to the speed of sound itself (something 
that will never be seen in everyday physics). 
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Chapter 2 

Models for heavy ion collisions 



2.1 Thermal model 



Since the transverse momentum spectra of partióles in heavy ion collisions are 
pretty similar to what is predicted by thermodynamical assumptions, they 
are sometimes described as coming from a source that is globally equilibrated 

El E [81 El [10]. 

The model starts with the phase space density, which is given as 

fUp) = ^ = ^^ = ^ l (21) 

^^"^'^^ dr d%d3p (27r)3exp(^)+a ' ^ ' 



Here, E = ^/p^ + m? is the energy, which will in the more general case of 
a moving source (moving with the four-velocity u^) be replaced by the scalar 
product p^u^. g and /i denote the degeneracy factor and chemical potential 
for the particle species considered, and T is the temperature. The latter is 
assumed to be the same for all species. a, finally, distinguishes the different 
spin statistics: For bosons (integral spins), it is «be = ~1) which results in 
the Bose-Einstein-distribution, whereas for fermions (half-integer spin) it is 
QíFD = +1, which leads to the Fermi-Dirac-distribution. For high energies 
and high temperatures both distributions become equal, a can be neglected, 
and one obtains the Maxwell-Boltzmann-distribution. 

A thermal model is by far the most macroscopic approach to heavy ion 
collisions. It can not account for local inhomogeneity and its applicability to a 
system exploding immediately and with high velocities is highly questionable. 
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2.2 Transport model 

An approach that accounts not only for local deviations from an assumed 
global symmetry, but really tries to model each particle and its trajectory 
through the space time, is the transport approach [Til [T2l [T3] . The evolution 
of the system foUows the relativistic transport (or Boltzmann-) equation 

p^9,/(x,p) = St{/} (2.2) 

Here, St {/} is called the coUision term, which contains Information on 
cross-sections and acts as source term for the density function. 

A set of Solutions to that formula is the equilibrium phase space density 
used for thermal and hydrodynamical models (12.11) . 

Transport models assume point like, classical particles whose mean free 
path is very large. The description of three-particle-collisions is very compli- 
cated. 

2.3 Lattice QCD 

An approach to solve the equations of QCD directly is to perform lattice 
calculations. Here, space and time are descretised (therefore "lattice"). The 
properties of infinitely vast matter in equilibrium, and therefore the Equation 
of State, can be studied [Hl [T51 [ISl [lU [TB]. It does not provide a dynamical 
description, but can give exact input to hydro- and thermodynamical models. 

2.4 Hydrodynamical model 

A model widely used to describe heavy-ion reactions is fluid- or hydrodynam- 
ics. It has been predicted as a key mechanism for the creation of hot and 
dense matter very early [Tül [2Ü] . Using this approach one assumes that the 
matter described is in local thermal equilibrium or at least in a state showing 
only small deviations from that. This assumption is not obviously met in a 
heavy ion reaction. Indeed, seen from any frame, the particle distribution 
functions of projectile and target are very different and far from being equili- 
brated. In the progress of the coUision, though, a locally equilibrated system 
may be formed. The part of the reaction before equilibration can therefore 
not be described by usual (1-fluid-) hydrodynamics. 
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One approach is to introduce several distinct fluids that each hold a equi- 
librated subsystem, e.g. one fluid for the target, one for the projectile and 
optionally a third for the evolving flreball. The whole system here is not 
equihbrated, but hydrodynamical evolution of each component is possible. 
Interaction between the components is implemented by considering the other 
components as source terms in the equations. 

This approach allows for modelhng of transparent nuclei, i.e. of the fact 
that at low energies the nuclei penétrate through each other. 

In so-called one-fluid hydrodynamics the beginning of the reaction can- 
not be modelled. Therefore, an additional model has to be apphed for the 
creation of the flrst equihbrated state, the so-called initial state. This can 
in principie be any kind of non-equilibrium model. The creation of the flrst 
equihbrated state can take different times, depending on the mechanism with 
which it is reached. Times of the order of Tq ~ 1 fm are reasonable. 

When the initial state is deflned, hydrodynamics start to work. Depend- 
ing on number and kind of the assumed symmetries in the initial state the 
fluid development may be solvable analytically or only numerically. Within 
this stage of the calculations, one has to assume an equation of state (EoS). 
It must also be chosen if the evolution describes a perfect fluid that is per- 
fectly equihbrated at any point or if small deviations are endorsed. In the 
latter case, several additional parameters like the heat conductivity, shear- 
and bulk-viscosity have to be introduced. 

In the case of a heavy ion collision the reaction zone is surrounded by 
vacuum, into which the matter will expand. Therefore, the energy-density 
and the particle density will decrease as a function of time. The system gets 
more and more dilute, and the assumption of local equilibrium at any point 
gets more and more unjustifled. A condition has to be deflned where the 
hydrodynamical evolution is stopped and something else is done. This step 
is called freeze-out. Whatever the condition is — speciflc lab- or proper- 
time, temperature, energy density, baryon number density — , it must be 
such that each part of the system crosses that condition and can be frozen 
out. The área in which freeze-out happens is usually a three dimensional 
hyper-surface, but it might have a flnite thickness, which would lead to a 
hyper-layer or hyper-volume. 

Apart from the cholee of the position of the surface, some other things 
may be adjusted for the freeze-out process. This starts with the kind of 
matter in the flnal state. Here, massive hadrons, massive quarks (without 
hadronisation so far), even massless quarks or hadrons may be assumed. 
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Quantum effects in the distribution function may or may not be taken into 
account. 

In any case, freeze out should conserve some quantities. Resides energy- 
momentum- and baryon number-conservation care has to be taken that en- 
tropy does not decrease. It is an additional constraint to require even the 
gross baryon number (number of baryons plus number of antibaryons) to be 
equal across the freeze-out process. 

Even for the surface it is not enough to state its position, its thickness has 
to be taken care of as weU. While the classical, standard approach considers 
it to be infinitely small, it may also be that freeze-out happens at a whole 
layer, which imphes that particles may leave the system that are cióse to, 
but not at the bordering condition. 

Starting with the Boltzmann transport equation (12.21) with a vanishing 
source term, one can define moments of the distribution as 

N^{x) = I ^P'f{x,p) (2.3) 

/i3 
-^p^p-f{x,p) (2.4) 

which are the baryon number current and the energy-momentum density, 
respectively. Since p^d^j,/ = O (see (12.21) ). it follows that 

a^iV^= O (2.5) 
d^T^'' = O . (2.6) 

These five equations (there is one for every u in (12.61) ) face fourteen inde- 
pendent variables whose time developments have to be found (four indepen- 
dent components of A^'^ and 10 independent components of the symmetric 
four by four-tensor T^'^). By assuming perfect local thermal equilibrium 8 
independent components of T^'^ can be eliminated. Then, one is left with 
6 parameters and only needs one more equation — the equation of state 
(EoS). This usually gives a relation between pressure, energy density and 
baryon number density. For a detailed description on how to decompose 
T^'^ into a ideal and non-ideal part and the different possibilities for various 
definitions please refer to [Hl [22] . 

Very simple equations of state consider ideal, ultra-relativistic (i.e. mass- 
less) gases. Here, the speed of sound is constant c| = 1/3. The pressure is 
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not dependent on the baryo-chemical potential. This EoS is considered good 
in the QGP-regime. For lower temperatures, the simplest case is a (massive) 
hadron-resonance-gas. Here, the speed of sound is given by c| = 0.15. 

The hydrodynamical model used in this thesis is the partióle in cell 
(PlC)-method which has been developed by Harlow and Amsden in the 
early 1960s [231 121] and upgraded to ultra-relativistic energies by Nix and 
Strottman in the 80s and 90s ¡251 ISSl [27] . It combines the advantages of 
fixed cells (Eulerian grid) and free moving particles (Lagrangian markers) 
and has e.g. exphcit baryon number conservation. An ideal gas with two 
flavours is assumed for the Equation of State. It has a constant speed of 
sound cs = 1/v^. The initial state has been developed recently by Magas, 
Csernai and Strottman and is exhaustively explained in [28]. Freeze-out hap- 
pens at constant time and goes from massless QGP to massive quarks with 
a restmass of rriq ^ 300 MeV. Hadronisation is abjured [29]. It conserves 
energy, baryon number and gross baryon number. 
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Chapter 3 

Jets and Médium 



In any coUisions of the kind 2 — í> 2 (i.e. two incoming particles produce two 
outgoing particles), the daughter particles will be, seen from the center of 
mass system, back-to-back-correlated. This is a simple consequence of con- 
servation of momentum. When the momentum transfer in a coUision 
becomes sufñciently large, the structures that play a role for scatterings be- 
come small enough not to resolve whole hadrons, but their constituents, the 
partons. In other words, at high energies partonic interactions domínate 
over hadronic interactions. In this región also perturbative QCD (pQCD) is 
applicable, because the strong coupling constant «s is small enough. 

When two partons scatter, their daughter particles usually fragment into 
hadrons after a short time, that means they créate new gg-pairs from the 
vacuum, thereby losing energy and form hadrons themselves as well as new 
hadrons, that typically go parallel to the leading hadrons. Theses bunches 
of hadrons are called jets. 

If nothing hinders the jets from being detected, as is the case in proton- 
proton coUisions, two jets will be seen that are usually exactly back-to-back- 
correlated. The center of mass of very high energetic partonic coUisions is 
usually the same as the center of mass of the proton-proton collision. 

When correlating the azimuthal angles (p of the observed hadrons one can 
recognize such events by two peaks at Acp = O and Aip = vr. This signal for 
di-jets is seen in pp-coUisions at high energies. 

In heavy ion coUisions, it is not so clear whether such a signal can be seen. 
The created jets might have to penétrate through the médium that is present 
in such reactions, except when the initial parton-parton collision happens at 
the surface of the médium and both daughter particles go tangential to the 
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Figure 3.1: Two particle correlations for p+p, d+Au and Au+Au-coUisions 
at i/J/vjv = 200 GeV. Trigger particles are in 4 GeV < p±Trig < 6 GeV, and 
associated particles are all those with a momentum bigger than 2 GeV and 
less than the trigger 2 GeV < p_Lassoc. < PxTrig- The curves show a clear 
suppression of secondary particles at Aip = vr in gold-gold-collisions when 
compared to the smaller systems. From ¡30]. 



surface. In all other cases, at least one of the jets has to go through the 
médium. Indeed, the most interesting case is a hard collision cióse to the 
surface with one daughter going right out of the collision zone and the other 
going the longest possible way through it. Here, studying the so-called near- 
side-jet (the one leaving the system immediately) can give insight on the 
parameters of the original hard collision and on the parameters of the away- 
side-jet. 

Depending on the properties of the médium the away-side-jet will look 
different to the experiment. For example, the cross-section for interactions 
between jet and médium might be very small, or the médium dilute enough, 
so that the jet goes through the médium almost undisturbed. This is the 
case at (comparatively) low energies. 

At central RHIC coUisions, however, there is no away-side-jet observed 
with an energy comparable to the near-side-jet [30], see figure [3A1 

Obviously, the away-side-jet loses its energy, which is in return absorbed 
in some way by the médium. Different models for this energy loss exist. 
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Figure 3.2: Two particle correlations for p+p, d+Au and Au+Au-collisions 
at y^SNN = 200 GeV with same threshold for trigger partióles as in figure 
13.11 but lower trigger for associated partióles 0.15 GeV < p^assoc. < 4 GeV. 
The data show sideward peaks at A(p ^ vr ± 1. From p7] . 



and the exact mechanisms are subject to current discussion in the field. A 
long time [211 E21 E31 EH E5] radiative energy loss has been considered the 
dominant mechanism. In 2003, Mustafa and Thoma [3E] re-considered energy 
loss by collisions, as had been predicted long before [371 EHl EHl HD] . Only 
recently, Mustafa and Thoma's results gained more attention jH]. Peshier 
¡l2l HSl un HS] has shown in 2006 that coUisional energy loss is indeed very 
important for the explanation of the magnitude of transport coefficients etc. 

In any case, secondary particles will have significantly lower momentum 
than the near-side-jet, which is a good explanation to the apparent complete 
disappearance of the away-side-jet in [30]. When lowering the threshold for 
secondary particles, one can see sideward peaks in the two-particle azimuthal 
correlations [16], see figure [321 Such a signal has been predicted as signature 
for mach shocks in 2004 [l8]. But also large-angle gluon radiation [l9l [50], 
jets defiected by radial fiow and cerenkov radiation [SU [52] are consistent 
with the observed away-side structure. 

Mach cones in nuclear matter have been predicted for cold nuclear matter 
[ini [531 [MI [551 [56], fermi liquids ¡53 M\ and QGP ¡591 [601 [61] and observed 
in heavy ion reactions at RHIC ¡62l[63l[6l[65l[66l[671[68l[69l[7Ü]. 

Recently, the predictions are supported by measured three-particle cor- 
relation spectra, where two azimuthal correlations are being opposed (see 
figure 13731) . Here, a clear distinction can be made between an indifferent en- 
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Figure 3.3: This sketch shows the angles used for three-particle correlations. 
Azimuthal differences from two different particles to the high-p^-particle are 
opposed. 



hancement of spectra by deflection etc. and mach cones or cerenkov radiation 
that pronounce exphcit directions. Peaks at (A0i, A02) = (vr ± 6, vr it fe) (on 
the bisector) would be present in all scenarios, they correspond to the fact 
that a jet typically consists of more than one particle, and a peak on the 
bisector merely says that two particles are emitted at the same angle. Peaks 
at (A01, A02) = (tt ± 6, TT =1= 6) (note the opposite sign!), on the other hand, 
show a correlation between markers on both sides of the backward direction. 

To distinguish between cerenkov radiation and mach cones, one can look 
at the p_L-dependence of the angle. In |51] it is shown that the angle of the 
cerenkov cone should be increasing very quickly with the momentum of the 
associated particles. 

Preliminary data from STAR (see figure 13.41) support the conical mecha- 
nisms (mach cones or cerenkov radiation), and deeper insight leave, according 
to [7T], little doubt on the dominance of mach cones over cerenkov radiation. 

Satarov [72] and Chaudhuri [73] have calculated that the jet angle in an 
expanding médium is also very dependent on the exact origin of the jet. In 
short, the angle rises when the jet does not come from the middle of the 
médium. 
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Figure 3.4: Three-partiele eorrelations from STAR, b (see text) is about 1.1 
radians = 63° here, whieh eorresponds (in statie médium) to a speed of sound 
of c| ~ 0.21 ± 0.8. See text for more details. From ¡71¡ . 

3.1 Calculated eorrelations in statie médium 

It is not trivial to calcúlate the expected angular distribution even for statie 
médium. The reason for that is that one has to use two different coordínate 
systems, both of which are spherical. One is the laboratory system, having 
the beam axis as the pole, the other one having the jet pole. Describ- 

ing the cone is very easy in the latter, which will be denoted {a, ¡3) for polar- 
and azimuthal angle, but the measured angle will be seen from the outside 
system {'&,íp). Both systems never are the same, because the jet can — in 
the model used — never go into the beam direction (see chapterU]). 

A jet propagating through statie médium will cause a cone appearing at 
a constant longitude, in the system where the jet is going towards the pole. 
The angular distribution is here 





or, in other words, the cone is at 




(3.2) 
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where the index a/3 defines the coordínate system to have the ¿-axis in the 
direction of the pole, which is the direction of the jet. Here, eos (vr — «o) = 
—eos («o) is fixed to the speed of sound (it is the usual mach angle), but 
/3 runs around the circle: /3 G [0; 27r). Each point on this circle is equally 
weighted with spectral enhancement, so that dN/df3 = const. Note that, as 
long as one stays consistent within one consideration, the exact choice on the 
¡3 = 0-direction is arbitrary, therefore, at the start, eos {¡3) and sin are 
interchangeable. A relative negative sign will occur if a jet in positive and 
negative beam direction are considered, because the direction of rotation of 
P has been changed. 

In order to obtain the spectrum dN/d(p as measured in a heavy ion ex- 
periment, one has to rotate the system a/3 to i!}ip. For a jet at mid-rapidity 
this is simple. In the following, the jet will always go in x-direction, where 
the azimuthal angle is zero (this is equivalent to using the difference angle 
A.(p). Then, the cone will be at 

(eos («o) \ 
sin (ao) eos (/5) . (3.3) 
sin («o) sin (/?) J 

Now, the azimuthal angle ip can be read off: 

y 

tan (üp) = — = tan («o) eos (/3) . (3.4) 

X 

dN/dip can be expressed as dN/df] ■ d^/dip, where the first factor is one (see 
above). The latter can be obtained from equation fl3.4p : 

/3 = — csinfi^l (3.5) 

(3.6) 



tan (ao) 
d(3 -(l + tan2((^)) 



d(p A/tan^ («o) - tan^ (ip) 

For the minus-sign in equation (13.61) refer to the statement above. 

In the case of a jet that is not going to mid-rapidity, the formulas become 
a lot more complicated. To characterize the direction of the jet, we take the 
angle r between the direction of the jet and mid-rapidity. Thus, r = 7r/2 — í9. 
The general idea is the same as before; we rotate fa/s (see equation (13.21) ) to 
f^^, solve for /5 and derive with respect to ip. 



3.1. CALCULATED CORRELATIONS IN STATIC MEDIUM 



21 



The rotation now happens with a more eomplieated matrix; we rotate 
around the y-axis with an angle of r: 

eos (r) O +sin (r) \ / eos («o) \ 
r^^ = { O 1 O sin («o) eos (3.7) 

—sin (r) O eos (r) J \ sin («o) sin (/5) y 

eos (r) eos («o) + sin («o) eos (/3) sin (r) \ 

sin («o) sin (/3) . (3.8) 

eos (r) sin («o) eos (/?) — eos (ao) sin (r) y 

As before, we read off tan (ip) as the ratio of y- and x-eomponent. Using 
Acos (0) + Bsin {(j)) = y/ + B'^cos {(j) — aretan {B/A)) and basie geometry 
one finds 



P = areeos 



— eos (r) tan {(p) 



tan («o) V ^ + (''") 



+ aretan ( — — — — - ) (3.9) 



sin (r) tan (if) 
and, after some more ealeulations 
d/3 1 + tan^ {(p) 



d(p 1 + tan^ ((/)) sin 



2 / \ {~sin (r) + eos (r) * 



tan^ («o) + tan^ (</)) [tan^ («o) sin^ (r) — eos^ (r 



.(3.10) 



Ahhough equation (13.101) seems very eompheated, two speeial eases ean 
be examined very easily. For r = O, whieh eorresponds to a jet in mid- 
rapidity, one ean re-obtain equation (13.61) . and for r = ±7r/2 one ean as well 
easily see that the result is =f1, whieh is expeeted beeause in these eases (the 
jet goes along the ¿-axis) ip and =F/3 are equal. The funetion is plotted for 
three different r, namely O, 30 and 60 degrees, in figure IXül 

For eertain angles r, the distribution fl3.10p is obviously undefined. In 
the eases where the root gets imaginary, there are no partieles emitted, the 
distribution should henee be set to zero. In the eases where the root diverges, 
a divergent measurement will be prevented by the faet that what is measured 
in an experiment is always the average of the distribution funetion over a 
small interval. 
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Figure 3.5: Calculated two-particle correlations in static médium. Shown 
are the correlations for a jet at mid-rapidity (solid line), for a jet at r = 30° 
(dashed line — 7r/6) and at r = 60° (dotted line — 7r/3) in the backward 
hemisphere. It can be seen how the máximum wanders "outside". 



However, at the border of the defined interval there is a peak. With the 
modulus of r, |r|, getting bigger, the defined interval gets bigger. For mid- 
rapidity it is P = (vr — arceos (cs) , tt + arceos (c^)), but for all other jets 
it will be larger. No jet will therefore contribute a peak between the mid- 
rapidity-peaks, but only outside. Therefore, the máximum of the distribution 
function will shift "outwards" , i.e. it will suggest a smaller speed of souncQ. 
Note that this effect always draws in one direction, and only a measurement 
at absolute mid-rapidity may reveal the true speed of sound — but then 
again, this all is for a static médium only. 



^Of coursc, this could be shown by integrating cquation (|3.10|) ovcr r and analyzing 
the resulting function. We rather take the figurativo approach. 



Chapter 4 

MACE — Mach Cones 
Evolution 

In order to study mach cones in a médium with realistic behaviour, a model 
has to be buih that either creates cones together with the médium (so that 
the cones become an inherent part of the evolution) or one that does add 
cones to the médium with hindsight, as perturbations. 

The first case requires the evolution algorithm to have a much higher 
spatial resolution than it would be reasonable to have for a hydro-code, else 
small (and localized) perturbations would be lost very quickly. Propagating 
a jet as a very strong perturbation within a locally equilibrated médium is at 
least a questionable thing to do. The arising dilemma would be eliminated 
if one only propagated the sound waves within the equilibrated system and 
considered the jet being an external source of energy and momentum that, 
if it is affected by the médium at all, only interacts outside of the hydro 
framework. 

In the latter case, on the other hand, backreaction towards the system 
has to be neglected. This is, though, a quite reasonable cholee for sound-like 
perturbations which are expected to be small (else they are not sound-like 
any more). But regrettably also energy-conservation is not fulfiUed. Neither 
the energy of the jet ñor the additional momentum which in the end will be 
calculated can be taken out of the system but have to be added. Therefore, 
one has to assume that the total energy is much bigger than the energy 
added. This, too, is not very hard to argüe for; even at RHIC energies 
{^/sKÑ = 200 GeV) adding a single 50 Ge V- Jet (which is tremendously high) 
would only change the total energy by 0.13 %. 
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Adding waves to the system after the evolution is calculated allows for a 
much higher resolution; the position of a wave can in principie be specified 
up to the precisión of floating-point operations at the calculating machine, 
whereas in the other case one was hmited to the grid size used for the system. 

Our approach is to take an existing hydro-evolution and impose a jet and 
the waves it creates after that as perturbations. 



4.1 Initialization and propagation of the jet 

Only the away-side-jet is considered, since the near-side-jet is not affected by 
and does not affect the médium. 

The jet starts from a random position within the médium at the first 
timestep calculated by the hydrodynamical code. This is, depending on the 
system considered, after few fm after the first collisions. 

The jet is created at a given point with a probability that is proportional 
to the energy-density at that point. The jet's direction is totally random. 
This is a reasonable thing to demand, since we consider an equilibrated sys- 
tem, in which collisions in any direction may occur. So, even with solid 
angle-dependent cross-sections one will get a spherical symmetric distribu- 
tion of secondary particle^. The only cut that is reasonably made is that 
we exelude all jets that will not end up in the detector. More specifically, 
we exelude for these studies all jets with a pseudorapidity \ri\ < 0.9, which 
is the acceptance of the ALICE-TPC [H]. This corresponds to an opening 
angle of A6 ^ 88.5° and to a covering of about 70 % of the solid angle. 

The jet is considered to be high energetic enough not to change its velocity 
and direction. Calculation of jet quenching is not within the scope of this 
thesis. It henee propagates in a straight line with speed of light through the 
médium and excites sound waves as it goes along. Its (final) direction will 
be used for correlation considerations in the end. When out of the médium, 
the jet does not excite sound waves anymore. 



^This does not apply to the final measured momentum distribution. but it must be 
true if we take each collision by itself. 
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4.2 The waves 

No premature assumptions on the shape of a resulting mach front can be 
made in an unforeseeable, inhomogeneous and non-statical médium. The 
reshaping of a mach front after boosting in and out of the fluid rest frame 
(FRF) and considering various ahgnments between the jet's and the fluid's 
direction have been discussed in several pubhcations [721 [75l [761 [73] , but none 
of them took into account a spatially inhomogeneous [721 [75] or a reahstic, 
three-dimensionally evolving [Tül [75] system. Also, they considered only the 
wave front. 

We take a different approach: we will propágate the single elementary 
waves and identify the wave fronts independent of the propagation. This 
allows for parts of the elementary waves to become a part of a wave front only 
after some time and possible deflection. Also, some parts may be swamped 
away out of the wave front. 

So, in order to model the sound waves we créate a lot of logical particles, 
so-called wave markers, at the position of the jet. The word "logical" refers 
to the fact that we do not assign any physical quantities to these markers 
yet; they only represent the position of the wave. The number of these 
wave markers per timestep is an adjustable (and numerical) parameter, in 
the standard setting it is nmarkers = 1 000. The entirety of the wave markers 
sent out at one timestep will be referred to as one elementary wave. The 
markers will be assigned random directions, so that after a short propagation 
in homogeneous médium one elementary wave should indeed be a spherical 
wave. Such an elementary wave is created at each timestep. Between the 
timesteps all "waves", i.e. all wave markers, are propagated. 

The propagation of the wave markers is straightforward: The only as- 
sumption made is that they move with the speed of sound relative to the 
fluid wherever they are. Their direction is adjusted by relativistically adding 
their initial velocity, v, to the flow-velocity of the underlying médium at the 
current point, Ü. Technically, this is done by adding the vectors correspond- 
ing to the respective rapidities and yu- 



V ® Ü 



(4.1) 





(4.2) 
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(normal vector- addition) (4-3) 

^ cxp(2|y/|)-l \ 
exp(2|y/|) + l \ 

vf. I í9/ = ^y, (4.4) 

When a marker crosses the border of a fluid cell, then its current propa- 
gation fl = f+ V ■ dt will be flnished with the oíd velocity, i.e. no deflection 
happens at the border of two cells. This is justifled because the timestep di 
is much smaller than the dimensions of the cells dx, dy and dz, so this will 
not happen all the time and will not cause too big an error. 

If a marker leaves the system, it is deleted. No particle emission at this 
point is assumed, and the wave is also not reflected. 



4.3 Freeze-out of the sound wave 

In the following we will discuss how to extract spectral enhancements from 
the position of the wave markers. This is not at all trivial. To explain 
why a certain method does or does not work, we will refer to the analytical 
test case of a static médium. Unless explicitly denoted, we claim that the 
argumentation is valid for non-static médium as well, but would be a lot 
more complicated to explain. 

4.3.1 Why trivial addition does not work 

A trivial way to evalúate the wave-markers is to assign a certain magnitude 
of perturbation to each of them and add this to the current cell's momentum 
and energy density. This may be done with a positive "delta" -function that 
adds a valué at the exact point of the marker, or with a smoother function 
that adds valúes at the position of the marker, but also before and after that; 
presumably, since we want to consider sound-like perturbations, we should 
add as much energy-momentum as we subtract at another place. The latter 
might also simúlate a kind of interference between different wave-markers. 

This way is a good physical cholee in order to obtain spatial Information 
about the perturbations. If such a spatial picture of the mach waves is to 
be made, this is the way to go. This is, in fact, what nature does, and what 
helps US to make photographs of aeroplanes or bullets (see flgureHTT]) showing 
a clearly visible mach cone. 
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Figure 4.1: A picture of a buUet moving through water. The buUet is moving 
faster than the speed of sound in the médium, therefore a mach cone develops. 
From ¡77¡ . 



Unfortunately, a picture of the mach cone is not asked for. We try to 
achieve a momentum distribution of the system, since this is what will be 
measured. The momentum distribution of a system altered using the above 
method will be equal to the unaltered distribution (or enhanced, if energy and 
momentum have only be added). Using constant time freeze-out (see section 
12.41) we intégrate over the whole "picture" and pick up the contributions 
along the mach cone that we saw on the picture. We could see it since their 
magnitude has been big due to the sheer number of wave markers in cióse 
proximity. But we also pick up the contributions opposite to it that do not 
contribute to the picture. They are not visible because there are only some 
at relatively big distances. Nevertheless, even if we would not discover their 
contribution on a spatial picture, their integral remains the same, no matter 
if they are cióse to each other or far apart, so they will cancel the contribution 
of the mach cone. 

Not even the interference one might build in will help here, because in- 
tegration is a linear operation. Taking all perturbations, squaring them and 
adding them to the system, which would only "count" the amplitudes and 
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resolve the "integral is constant"-problem, will lose all information about the 
waves direction. These effects are not unphysical and exactly what would 
happen with the air perturbed by a supersonic plañe. The problem arises 
with the heavy-ion-way of measuring. When studying a plañe we take a 
picture and look at where the mach cone is. This means we look at the 
amplitude. If we do not take several pictures, we have, in fact, no idea on 
where the sound wave will be at a later time step. Another way on how to 
watch a fast plañe is to stand still and listen. We can hear the perturbation. 
The analogous thing to do in the hydro-evolution is to go to one place and 
wait for the wave to come — i.e. freeze out at constant place. Due to the 
assumption that particles will not change their direction after freeze-out spa- 
tial information would directly be mapped to momentum information. This 
method fails for the used freeze-out method (see again section [2^ page [T3!) . 

4.3.2 A coalescence-like model 

Following the idea that in the wave front many wave markers are cióse to 
another and going into the same direction, it might be a good idea to only 
consider markers that satisfy this condition. It would require two new param- 
eters; the distance in which a different marker is looked for, and a maximal 
angle by which the directions of the two markers may deviate. 

Unfortunately, one finds that in the analytic test-case the distance be- 
tween two markers emitted into the same direction a at to subsequent time- 
steps, will be dr = dt -^/l + c| — 2c5Cos(q;), which will be minimal at a = 0. 
Forward markers can therefore not be excluded by this method, although they 
do not contribute to the mach cone. 

This ansatz anyways has not much to do with interference; markers in 
short distances in the direction of their movement will rather interfere than 
contribute to collective movement. Markers that are next to each other, 
though, will contribute. This fact might be taken into account by not using 
the distance of two markers as criterion. Instead, one can model the distance 
in a different way: 

Usually, one can reinterpret the concept of distance between r and rV as 
taking a sphere around the point r and adjusting the radius so that r/ is on 
the surface of that shell. The radius dgph is then the usual — spherical — 
distance between the points. 

The idea that is to be presented here is not to take a sphere and adjust 
its radius, but to take an ellipsoid with given excentricity and adjust the 
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remaining parameter (one of the semi axes) so that the second point is on 
the shell of the elhpsoid. 

For these assumptions it is reasonable to have the semi-axes of the elhp- 
soid ahgned along the direction v oí the considered wave marker; the elhp- 
soidal should be symmetric in the azimuthal direction with respect to v. In 
fact, the semi-axis paraUel to v is, for the reasons stated above, the shortest. 

The (imphcit) equation for an eUipsoid satisfying the above constraints 
having semi-axes A orthogonal and B longitudinal to v is 



the spherical distance for B ^ A, i.e. lim^^i d^w = ásph- 

In the same case as above we can now calcúlate ddi as a function of a. 
It turns out that we indeed have a máximum for cícii at forward directions 
now instead of a mínimum, and a mínimum at eos (a) = ±i?^ cs/{B'^ — A^), 
which is, though, always smaller than the expected angle eos (a) = cs- To 
reach it, one would have to set A = O, which is not reasonable. Using this 
ansatz would henee always underestimate the angle. 

4.3.3 What else can be done: The Mantle Method 

There is a way of looking at a picture of a mach cone and predicting where 
it will be in the next instant; we consider the movement to be perpendicular 
to the wave-front. Finding out the directions perpendicular (or parallel) to 
the front is easy to do by hand, one "sees" where the front goes. It is not so 
easy to find an algorithm that discovers these directions automatically. 

Our approach is to try and reshape the región occupied by wave markers 
with a set of lines, referred to as wave-lines. This surface will be interpreted 
as the wave front. 

The wave lines are equally distributed azimuthally around the jet's di- 
rection V. Their nodes are positions of marker particles. The rough idea is 
to go in opposite direction of the jet, starting at its position, so that seen 
from the jet it looks like standing at the pole of the earth and considering 
the longitudes going away covering the surface of the earth (apart from the 
earth being a sphere, not a cone). 




(4.5) 
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For the construction of each wave line, a subset of all wave markers is 
considered. A slice azimuthal to the jet axis is taken into account for every 
wave line. The center of the azimuth, though, is not necessarily on the 
trajectory of the jet, but is the geometrical center ("center of mass", though 
no masses are present) of a given elementary wave, so it changes for markers 
from different elementary waves. The azimuthal slice is classified by an angle 
4>iine- It is the angle between the plañe spanned by v and the wave line on the 
one side and the plañe spanned by v and á on the other side, where á is some 
arbitrary, fixed direction perpendicular to v. We shall use the convention 

O = 77 TV (4.6) 

\v X z\ 

b = vxá , (4.7) 

where z points in the z-directior0. We cali the right-handed coordínate sys- 
tem V, á, b the í;a&-system. 

Now, we can classify the wave markers by assigning them a similarly 
defined angle, </3marker- In principie, this angle is no different from (f>\inc, but we 
allow the intersection line between the two respective planes to change: While 
the v-á plañe remains a constant reference, the second plañe is now spanned 
by V and the difference vector of the wave marker fj and the geometrical 
middle tqu of its elementary wave rdm = rj — rcM (as indicated above). It 
should be noted that for these considerations it is important to track the sign 
of the angle. 

This problem can be composed into a projection of f¿ifí onto the á/b- 
plane and the angle between that projection and á. Considering rdifj = 
Aá + Bb + Vv, the projection is r± = Aá + Bb. The angle is then 

V'marker = arCCOS |^--^=^=^ SgU (B) (4.8) 

where A = {fj — rcivi) ■ o and B = {fj — rcivi) ■ b are the projections of rj_ 
towards á and b. 

For a given wave line with we define a new coordínate system vfg((j)) 



^This choice for a is practical since wc do not considcr cvcnts with the jet going into 
the longitudinal direction z, so we do not have to check whether á is non-zero. 
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which is defined by 

eos (0)iine a + sin b and (4.9) 
vx f . (4.10) 

Furthermore, all wave markers with v^^arkcr ^ [0imc — ^0/2; 0iinc + ^0/2) 
are considered as nodes, where A0 = 2n/Ni[^cs is the angle between two 
wave lines. 

The sheer surface of the área might be very distorted or irregular due to 
numerical reasons. For instance, due to the finite number of timesteps, the 
surface will contain large portions of the same elementary waves (in the hmit 
of infinitely many elementary waves each of these only contribute one point, 
which is a circle in a plañe orthogonal to the jet's direction). In order to 
account for this numerical problem, one has to maximize the angle between 
the wave line and the jet axis when looking for the next node. This angle, 
denoted a, between a node and a candidate fj has to be maximized. It 
can be written analogous to (14.81) : 

a =^ arceos ( , ^ | sgn (F) (4.11) 

where —V and F are the projections of — fj onto —-O and /, respectively. 

Another numerical aspect is the random direction of the wave markers 
and their finite number. It might be the case that in the área where a given 
elementary wave should represent the surface there is no wave marker to 
determine a possible node. Therefore, one has to give the possibility to skip 
a single elementary wave and go straight onto the second next one without 
any node taken from the next. AUowing for more than one elementary wave 
to be ignored may conceal physical effects when collective fiow is going into 
different directions at different points. Figure 14.21 sketches the mentioned 
cases and shows the surfaces without and with the named improvements. 

The search for the next node holds one additional pitfall: since the pos- 
sible wave markers are not all in one plañe, a slight deviation might occur 
and the wave line might finally drift out of its domain; in the worst case, all 
lines would end up with the same nodes after some distance. Therefore we 
have to ascertain the azimuthal position of the wave lines. Henee, the nodes 
do not correspond to actual positions of wave markers, but to their rotation 
onto the fixed ¿//-plañe. Since we can assume that the number of lines is big 
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Figure 4.2: A schematic view on how the Unes are constructed, assuming that 
all wave markers are in a common plañe. In (a) the basic smoothing is shown 
(maximizing the angle). In (b) the effect of leaving out single timesteps is 
shown. Finally, in (c) a dynamically created wave line from data in static 
médium is shown. 

The dashed line shows the jet-trajectory, the dotted ares the elementary 
sound waves, the dash-dotted lines the wave lines how they would be without 
optimization and the solid lines the final wave lines. 

enough to provide small angle approximations for d0, we do not rotate, but 
merely project the vector onto that plañe. From a candidate fj we therefore 
get a new node rj+i with 



When looking for the next node, we have to exelude fj. To make the algo- 
rithm faster one can exelude all wave markers that lie "behind" the current 
node (that means, closer to the position of the jet). This can be checked by 
the relation {fj — fi){—v) > O, from which foUows 




(4.12) 



TjV < TiV 



(4.13) 



When all nodes are collected, particle fiow is considered to be orthogonal 
to the wave lines. The lines are taken as is. To find the places on which 
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Figure 4.3: Particle flow is assumed to be at fixed distances (the length of 
the wave line is measured) perpendicular to the current connection between 
next and previous position. No additional smoothing is done. Also it is not 
compensated if two or more nodes lie between two "freeze-out-points" . 



to insert momentum we go along the wave lines, starting from the jet, and 
insert momentum after fixed distances. Momentum is inserted in the v/f- 
plane orthogonal to the wave line. No special treatment for momentum 
insertion at a node is done; orthogonality refers here to the part of the wave 
line before the respective node (see figure |4]3¡). 



4.4 Energy-conservation 



As indicated before in the beginning of this chapter (page[23l). energy cannot 
be conserved in this approach. When counting the positions of momentum 
insertion before defining the amplitude of insertion one can, after measuring 
the distance l the jet has been going already, set the total amount of energy 
inserted AE to 



da; 



(4.14) 



where dE/dx\^^^ is taken from a different calculation. This will assure that 
the energy deposited in the mach región is monotonous increasing with time 
and has reasonable valúes. 



34 



CHAPTER 4. MACE — MACH CONES EVOLUTION 



4.5 Perfect Background Subtraction 

When the system is altered to contain the information of the sound waves, 
the momentum space distribution (and in principie the whole phase space 
distribution) can be evaluated to give interesting valúes. Actually, for two- 
and three-particle correlations, background from radial, directed and elliptic 
flow has to be subtracted in order to see a signal. For experiments, back- 
ground determination and subtraction is a major task. 

In the theoretical approach taken, it is possible to do two different things: 
One could take the original system and try to reproduce the analysis that has 
to be done in experiments. This would allow for better direct comparison to 
experimental results, but the data would be specific to a special experiment 
and background subtraction method. Comparison to data extracted with 
future, possibly better, methods would not be very easy. 

The other way is what we cali "perfect background subtraction". The 
idea is to leave out the background from the beginning. In this method, 
the system is not, as indicated in section 14.3.31 altered, but the information 
about where the system would be altered — the clean signal — is evaluated 
directly. This will result in the cleanest possible signal which would represent 
the limit of experimental evaluations. 

Here, we also have no need for correct estimation of jet energy loss, since 
we only need to consider that all alterations of the system have the same 
magnitude. 



Chapter 5 

Results from MACE 



There will basically be four different kinds of plots used in this chapter, which 
are now explained briefly. 

First of all, two-particle azimuthal correlations dN/d{Aip) vs. d{A(f) will 
be presented, where the near-side jet usually creates a máximum at Aip = 0. 
This peak will not be visible in our plots, since in MACE there is no near-side 
jet. The graphs in those plots show how many partióles have been present 
at an angle ipi = ip^ct + Aip. Since I have no particles in my model, the y- 
axis here is in arbitrary units, it counts insertions in the respective direction. 
Cmp. figures IXT] and IX^ 

A little more detailed insight to the actual situation might be taken in a 
three-dimensional two-particle correlation d'^N/d{Aip)/d{A-&) vs. d{Aif) and 
d(A'í9). On the y-axis, there is a second independent variable, the difference 
polar angle Ai? = -i^j — í^jct, analogous to d{A(p), which is again on the x-axis. 
The colours indícate how many particles have been at the respective angles 
with respect to the jet. 

Very similar is the three-dimensionally plotted two-particle correlation 
d"^ N / d{Aip) / d{Ari) vs. d{Aip) and d(A?7), where d(A?7) is the difference in 
pseudorapidity between jet and particle. 

Both these three-dimensional two-particle correlation only give new in- 
sight for an event-by-event analysis. If data are accumulated over many 
events and then correlated, they become no more meaningful as the (easier 
to understand and produce) two-dimensional plots. 

More insight, even over a sample of many events, can be taken with three- 
particle correlations dN"^ / d{Aipi)d{Aip2) vs. d{Aipi) and d{Aip2)- Here, the 
colours ( "¿-axis" ) denote how many particles have been present at an angle 
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m = 100, 1 = 72, d = 2 





m = 10000, I = 72, d = 2 



012345 012345 

A9 A<p 



Figure 5.1: The influence of the number of wave markers m on two-particle 
correlations (in static médium, after 10 time steps). From left to right, the 
number of wave markers goes up by a factor of 10 between the graphs. While 
the raising of m from 100 (left) to 1000 (middle) gives a much better signal, 
the improvement between 1000 and 10 000 does not change the quahty a lot. 



of Aifi while another partióle has been present at A(f2- Given that the 
correlations are done on an event-by-event-basis, a signal here will not be 
concealed or artificially created by adding the data from many events. Cmp. 
figures 13.31 and 13. 4[ 

The impact parameter is always given in fractions of the máximum pos- 
sible impact parameter, which is the sum of the radii of the involved nuclei. 
E.g., for gold-gold coUisions, 6 = 10 % = 0.1 ■ (2 ■ Rau) ~ 1-36 fm. 



5.1 Characterizing the algorithm 

The first thing that should be done with a new model or algorithm is to 
test it with an analytic test case or a known approximation (whichever is 
applicable), and, in the case of an algorithm, to test it against its behaviour 
when changing unphysical, numerical parameters. 

The analytic test case I am referring to is a static médium, i.e. a system 
that has a coUective fiow velocity u = O everywhere and at all timesteps. 
The algorithm is the one described in chapter IH especially in sections 14. H 
OI gXSl and 

When taking the approach as described in section l4l5l namely not adding 
momentum to a hydro system and freezing out afterwards, there are three 
numerical parameters left in the MACE-model: The number of wave markers 
created at each timestep m = rimarker, the number of wave lines / = nunes and 
the number of momentum insertions d = ripoints per one di length of a wave- 
line, where di is the length of one timestep. 
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m = 1000, l = 36.d = 2 




m = 1000, 1 = 72 

; 


, d = 2 


m = 1000,1 = 360 


,d = 2 
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Figure 5.2: The influence of the number of wave Unes l on two-particle corre- 
lations (in static médium, after 10 time steps). From left to right, the number 
of wave Unes goes up. While the raising of / from 36 (left, corresponds to 
an angle of A^ünc = 10°) to 72 (middle, A^unc = 5°) gives a much better 
signal, the improvement between 72 and 360 (A^ünc = 1°) does not change 
the quahty a lot. 



With static médium and no thermal smearing being done, the result- 
ing signal with perfect parameters should be very sharp. With imperfectly 
shaped surface, the angular distributions will be much wider, since, when 
the lines are jagged, the normáis will be in arbitrary directions. In the two- 
particle correlations shown in figures 15.11 15.21 and 15. 3^ the expectation is not 
to get two delta-functions. Instead, between the peaks the correlation func- 
tion should behave similar to an inverse circle, see section 13.11 and equation 
(13.61) . Outside of this región, it should indeed be zero. 

Figures 15.11 15.21 and 15.31 show two-particle correlations dA^/d(A(/9) for 
different valúes of m, l and d. The plot in the middle is the same in all 
figures; it is made with the valúes that are actually used. In 15.11 and 15.21 it 
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Figure 5.3: The influence of the density of momentum insertions d on two- 
particle correlations (in static médium, after 10 time steps). From left to 
right, the density of momentum insertions goes up. Neither step seems to 
change the plots a lot. 
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can be seen that the quality improves a lot from the worst chosen resolution 
to the middle one, but only minor changes appear above that. The signal 
seems to be rather independent on the density of momentum insertions d 
(cmp. figure [573]) . 

It should be emphasized what the parameters change: A bigger m makes 
the wave hnes smoother, while / and d increase the total number of momen- 
tum insertions N. Correlation is an (A^^)-algorithm, so both / and d should 
be as small as is reasonable. 

All in all, it can be concluded that the algorithm shows convergent be- 
haviour with respect to the parameters, so finite valúes can be taken. In the 
foUowing subsections, unless differently denoted, all evaluations are made 
using the valúes m = 1 000, / = 72 and d = 2. 

5.2 Mach cones at RHIC 

In the foUowing section, we discuss the results for correlation functions in 
gold on gold-coUisions at ^/sÑn = 130 GeV for different impact parameters. 
The influence of jet triggers (both direction and origin) is examined in section 
15.31 on central LHC-collisions. 

5.2.1 Central collisions 

The central RHIC data, when averaged over 1 000 different jet origins and 
directions, yield a double-peaked two-particle correlation function and two 
off-diagonal peaks in the three-particle correlations (see figure 15.41) . The 
máximum of the correlation, however, is shifted from the mach angle of 
a ~ 0.96 "outwards" by about 6a ~ 0.2 radians. The corresponding speed 
of sound to this angle is c^' ~ 0.17. This valué is very cióse to the 
expected result for a hadron-resonance gas of Cs ^0.15. 

5.2.2 Mid-central collisions 

As the central RHIC data, also the mid-central RHIC data (for impact 
parameter 6 = 25 %) show two (four) peaks in the two-(three-)particle- 
correlation function (see figure [575]) . Here, the same deviation from the per- 
fect mach angle is found, which is again much closer to the mach angle of a 
hadron-resonance gas. 
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Figure 5.4: Two- and three-particle correlations obtained for central RHIC 
coUisions (6 = 0%, red solid line (above) and solid lines (below), and b = 
10 %, blue dashed line (above) and crosses (below)) at ^/sÑÑ = 130 GeV. 
The dependent variable is given in arbitrary units. 



40 CHAPTER 5. RESULTS FROM MACE 



Au+Au 65+65 AGeV mid-central 





\ V 


r L 

I / \ 


rri 










\ \ \ 
















/ 


\ / \ 


■D 








1 2 


\ \ \ 

3 4 5 


5 












4 






3 






9- 












1 









-dN2/d(A(pi) d(A(p2), a.u. 

lili 



o 1 2 A(pi 3 4 5 



Figure 5.5: Two- and three-particle correlations obtained for central RHIC 
collisions {h = 25 %) at y/sÑÑ = 130 GeV. The dependent variable is given 
in arbitrary units. 
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Figure 5.6: Two- (above) and three- (below) particle correlations for central 
lead-lead-collisions at a/saT/v = 5.5 TeV. These figures show the sum over all 
correlations from 1 000 events with arbitrary jet orientation and origin. 



5.3 Mach cones at LHC 

In the foUowing section, we discuss the results for correlation functions in 
lead on lead-collisions at y/s^N = 5.5 TeV for central coUisions. 

Correlation functions measured for LHC energies resemble the picture 
from RHIC. Again, the "minimum jet-bias"-data show a shifted máximum 
with respect to the mach angle. It is also at A0max. ~ 1-15 rad, which 
corresponds to a speed of sound of c| ~ 0.17 ± 0.07. The minimum of 
the correlation function is about 30 % of the maximal valué, as has been 
before. Figure 15.61 (lower part) proves the existence of two máxima in most 
events (the two off-diagonal peaks). In figure 15771 two-particle correlations 
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Figure 5.7: Two-particle correlations for central lead-lead-coUisions at 
y/sÑÑ = 5.5 TeV as function of Aip and A-d (a) and A-q (b). These fig- 
ures show the same data as figure 15. 6[ It can be easily seen that these plots 
do not show the expected ring due to the average over so many events. 
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O 1 2 Acp 3 4 5 

Figure 5.8: Two-particle correlations for central lead-lead-collisions at 
^/sÑÑ = 5.5 TeV for two different classes of jets: The red solid line shows 
the result for 0.4 < ?7jct < 0.9, the blue dashed hne for |?7jet| < 0.1. 



are plotted for the same set of events, with spatial resolution now also in 
polar direction. This figure shows why this way of plotting a correlation 
function is not very helpful; there is no ring structure visible if one considers 
many events at the same time. All Information that can be obtained from 
these plots is that there is actually never a signal near the near-side-jet, and 
finite correlations at Aip O come from far away polar positions. The rest 
is contained in the two-dimensional two-particle correlations already. 

In an experiment, several triggers can be applied to the jet. While a 
trigger on the azimuthal position of the jet does not make sense in the az- 
imuthally symmetric central collisions considered here, it is interesting to take 
a look at different rapidities of the jet. In figure [5^ the two-particle correla- 
tions for forward jets (0.4 < ?7jet < 0.9) and mid-rapidity jets (jr/jetl < 0.1) 
are shown. This is a case where in static médium there would be a visible 
difference. Still, in the moving matter, both correlations look rather similar. 
Even the mid-rapidity data suggest a different speed of sound than is actually 
present. The error bars are the same as before. 

A much harder task on the experimental side is to find out where a jet 
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Figure 5.9: Two-particle correlations for mid-rapidity jets in central lead-on- 
lead-coUisions at ^/sÑÑ =5.5 TeV. 



carne from. Although this is not directly measurable, the effect of different 
origins can be studied with MACE. The two-particle-correlation functions for 
three jets going in mid-rapidity that started next to each other in the center 
of the colhsion and left as well as right to that (at 70 % of the way out of the 
médium) are opposed in figure 15.91 The centrally started jet, finally, shows 
the behaviour expected in static médium, a symmetric correlation function 
(symmetric around Aif = vr, that is), that peaks at A(p ^ tt ± 1. It is 
also revealed on this figure how much the starting point of the jet affects 
the outcome of the mach cone. The left and right jet are, though, pretty 
mirror-symmetric to each other, having one peak at the mach position and 
the other at about Aip ^ ±0.7. 

Concluding the results, one can say that the double-humped structure is 
a very robust signal. If there are mach cones in LHC coUisions, they should 
be visible. The máxima, though, are not at the naively expected angle. Only 
mid-rapidity jets with symmetric correlation functions can give insight on 
the speed of sound present when the cone was formed. In all other cases, 
the speed of sound will be underestimated. Due to the steep behaviour of 
the cosine function at the interesting valúes around 1 radians, even small 
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deviations from the perfect mach angle will give completely different results 
for the speed of sound. Care should therefore be taken to read off a speed of 
sound from a measured correlation function. There be dragons. 
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Chapter 6 

Summary and Outlook 



The properties of hot and dense nuclear matter have been subject to a lot of 
research projects in the last decades. To study this kind of matter, atomic 
nuclei (so-called heavy ions) are coUided at very high center of mass energies. 
From the observable particle spectra researchers try to get information about 
the different properties of matter. 

After the publication of some experimental results of the currently most 
powerful heavy ion-accelerator, RHIC, it became probable that hydrodynam- 
ical descriptions of nuclear matter in such collisions work pretty well. 

In Hydrodynamics, the spread of sound waves is intimately connected to 
the Equation of State (EoS). To get information about this it is good to find 
a signal for (or from) sound waves. 

A high energetic parton travelling through the médium {a. jet) might loóse 
its energy by bouncing off other partons. Its energy will then propágate 
through the médium as sound waves. With an ultra-relativistic jet, i.e. one 
moving practically with the speed of light, those sound waves would — since 
the speed of sound is smaller than the speed of light — interfere to form a 
cone, the so-called mach cone, as is the case with an airplane with supersonic 
velocities. 

In the presented thesis the model MACE (MAch Cones Evolution) is de- 
veloped, which simulates the propagation of sound waves through a médium 
that itself moves with relativistic speeds in the order of magnitude of the 
speed of sound. 

In Chapter H] it is shown how this model can be used to simúlate the 
sound waves of a moving source and automatically find collective effects as a 
mach cone ("supersonic boom") and analyse them. Here, the input needed 
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is the velocity-distribution of a médium as a function of space and time, in 
order to propágate the waves with a constant speed with respect to the local 
rest frame of the médium. 

The elementary spherical waves are modelled by logical points. Their 
entirety represents the sound wave. To recognize the mach cone automati- 
cally, the space región in which these points are, is mantled, i.e. the surface 
is determined. That means that interference-effects are not modelled, but a 
geometrical procedure is used. The surface is converted into a signal directly, 
it is abdicated to first créate an altered speed profile, as it is accessible to 
the experiments, to re-obtain the signal from that. 

From normáis to this surface a spectrum is calculated, which is measured 
as increased particle emission in the experiments. 

MACE is tested for its behaviour when changing numerical parameters in 
section lSTTl In sections [5 .21 and [5731 MACE is applied to the speed profiles from 
gold-gold-coUisions at y/sÑÑ = 130 GeV (RHIC) and lead-lead-coUisions at 
■y/sÑÑ = 5.5 TeV (LHC) and the results are presented. 

It turns out that, when averaging over many jets, no significant differences 
between the correlations from different starting points can be seen. Indeed, 
the added-up correlation functions are indistinguishable from those expected 
from static médium. Even when considering different special cases on the 
direction of the jet, the picture does not change. Only when triggering on 
the origin of the jet, which is experimentally the hardest of all tasks, one can 
see a pronounced change in the shape of the correlations. Still, for centrally 
started jets, there is no qualitative difference between LHC-data and static 
médium. 

In all averaged cases considered, the apparent, measured speed of sound is 
smaller than what has been input to the calculation and what has been used 
to propágate sound. This effect is strong enough to simúlate a hadron gas 
with c| = 0.15 in favour of the (actually present!) quark-gluon-plasma with 
c| = 1/3. As stated above, only mid-rapidity jets starting from the middle 
of the médium (recognizable by rather symmetric correlation functions) can 
reveal the true speed of sound. It is emphasized that an attempt to read of 
the speed of sound from an (in terms of the trigger) unqualified correlation 
function is doomed to fail. 

Since the underlying hydrodynamical algorithm (the PIC-Code, see sec- 
tion 12. 4p currently only uses an EoS with massless Quark-Gluon-Plasma, it 
could not be studied how the data change with variable speed of sound. The 
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porting of MACE to hydro-codes that do use a more complicated EoS and the 
changing of the speed of sound depending on the properties of the underlying 
médium will be a future project that may give interesting results. 

The next point that may be done with MACE is comparing the results 
from this approach to the data obtained by inserting a jet during the creation 
of the hydro-evolution. Here, differences between those models could be 
studied in detaiL 
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C.l Deutsche Zusammenfassung der Diplo- 
marbeit — Germán Abstract to the the- 
sis 

Die Eigenschaften von heiBer, dichter Kernmaterie sind Gegenstand vieler 
Forschungsvorhaben der letzten Jahrzehnte. Um solche Materie zu unter- 
suchen, werden Atomkerne (sog. Schwerionen) bei extrem hohen Schwer- 
punktsenergien zur KoUision gebracht. Aus den dabei messbaren Teilchen- 
spektren versuchen die Forscher, Aufschlüsse über verschiedenste Eigenschaf- 
ten der Materie zu erlangen. 

Nach Veroffentlichung einiger Messergebnisse des momentan leistungsfáh- 
igsten Beschleunigers für schwere Atomkerne, RHIC, scheint es wahrschein- 
lich, dass hydrodynamische Beschreibungen für Kernmaterie in solchen Kol- 
lisionen recht gut sind. 

In der Hydrodynamik ist die Ausbreitung von Schall eng mit der sog. 
Zustandsgleichung der Materie verbunden. Um Aufschluss über diese zu 
erhalten, ist es daher gut, ein Signal für bzw. von Schallwellen zu finden. 

Ein sich durch das Médium bewegendes hochenergetisches Parton (ein 
Jet) kann eventuell durch StoBe mit anderen Partonen Energie verheren, die 
dann in Form von Schallwellen durch das Médium propagiert wird. Bei 
einem ultrarelativistischen Jet, also einem, der sich praktisch mit Licht- 
geschwindigkeit bewegt, würden sich — da die Schallgeschwindigkeit kleiner 
ais die Lichtgeschwindigkeit ist — solche Schallwellen zu Kegeln (sog. Mach- 
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kegel) überlagern, wie bei Flugzeugen mit Geschwindigkeiten jenseits der 
Schallmauer. 

In der vorliegenden Arbeit wird das Modell MACE (MAch Cones Evolu- 
tion — Mach-Kegel-Entwicklung) entwickelt, das die Propagation von Schall- 
wellen durch ein Médium simuliert, das sich selbst mit relativistischen Ge- 
schwindigkeiten in der GroBenordnung der Schallgeschwindigkeit bewegt. 

In Kapitel H] wird gezeigt, wie dieses Modell benutzt werden kann, um 
die Schallwellen einer sich bewegenden Quelle zu simulieren und kollektive 
Effekte wie einen Machkegel ( "Uberschallkegel" ) automatisch zu finden und 
zu analysieren. Hierbei wird ais Input die Geschwindigkeitsverteilung eines 
Médiums ais Funktion von Ort und Zeit benutzt, um die Wellen mit kon- 
stanter Geschwindigkeit im Ruhesystem des Médiums zu propagieren. 

Die elementaren Kugelwellen werden dabei durch logische Punkte mod- 
elliert, deren Gesamtheit die Schallwelle darstellen. Um den Machkegel au- 
tomatisch zu erkennen, wird die Raumregion, in der sich diese Punkte aufhal- 
ten, ummantelt, d.h. die Oberfláche wird bestimmt. Das bedeutet, dass nicht 
Interferenz-Effekte modelliert werden, sondern geometrische Uberlegungen 
benutzt werden. Die Oberfláche wird direkt in ein Spektrum umgewandelt, 
es wird darauf verzichtet, zunáchst ein verándertes Geschwindigkeitsprofll zu 
erstellen, wie es den Experimenten direkt zugánglich ist, um dann das Signal 
daraus zu ermitteln. 

Aus Normalen auf dieser Oberfláche wird dann ein Spektrum ermittelt, 
das im Experiment ais erhohte Teilchen-Emission gemessen wird. 

MACE wird in Kapitel [5] zunáchst auf sein Verhalten gegenüber der 
Veránderung von numerischen Parameter getestet (Abschnitt 15.11) . In den 
Abschnitten 15.21 und 15.31 wird MACE dann auf die Geschwindigkeitsproflle 
von Gold-Gold-KoUisionen bei ^/s^ = 130 GeV (RHIC) und Blei-Blei- 
Kollisionen bei a/s/vtv = 5.5 TeV (LHC) angewandt und Ergebnisse werden 
vorgestellt. 

Es zeigt sich, dass, wenn über genug Jets gemittelt wird, keine signiflkan- 
ten Unterschiede zwischen den Korrelationen von verschiedenen Startbedin- 
gungen sichtbar sind. In der Tat sind die (aufaddierten) Korrelationsfunktio- 
nen nicht zu unterscheiden von denen, die bei statischem Médium erwartet 
werden. Sogar, wenn man verschiedene Spezialfálle für die Richtung des Jets 
betrachtet, ándert sich das Bild nicht. Nur bei triggern auf den Ursprung 
des Jets, was eine sehr schwierige experimentelle Aufgabe ist, kann man eine 
ausgeprágte Veránderung in der Form der Korrelationen sehen. Doch wieder 
ergibt sich für zentral gestartete Jets keine qualitative Veránderung zwischen 
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LHC-Daten und statischem Médium. 

In alien betrachteten, gemittelten Fallen verschieben sich die Korrela- 
tionsfunktionen so, dass man naiv eine Schallgeschwindigkeit messen würde, 
die um etwa das eineinhalbfache kleiner ist ais die tatsáclilich benutzte Schall- 
geschwindigkeit. Dadurch kann die Existens einer anderen Form von Materie 
(Hadronengas statt Quark-Gluonen-Plasma) vorgetáuscht werden. Wie oben 
angedeutet, konnen nur Jets in Mittrapiditát, die in der Mitte des Médiums 
starten (diese Jets sind erkennbar an einer symmetrischen Korrelationsfunk- 
tion), die wahre Schallgeschwindigkeit aufzeigen. Es wird unterstrichen, dass 
jeglicher Versuch, eine Schallgeschwindigkeit an einer Korrelationsfunktion 
abzulesen, zum Scheitern verurteilt ist, wenn man nicht genau weiB, auf was 
getriggert wurde oder der Trigger nicht auf das Problem passte. 
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